Modeling the Relaxation of Polymer Glasses under Shear and Elongational Loads 
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Glassy polymers show "strain hardening" : at constant extensional load, their flow first accelerates, 
then arrests. Recent experiments under such loading have found this to be accompanied by a striking 
dip in the segmental relaxation time. This can be explained by a minimal nonfactorable model 
combining flow-induced melting of a glass with the buildup of stress carried by strained polymers. 
Within this model, liquefaction of segmental motion permits strong flow that creates polymer-borne 
stress, slowing the deformation enough for the segmental (or solvent) modes to then re-vitrify. Here 
we present new results for the corresponding behavior under step-stress shear loading, to which 
very similar physics applies. To explain the unloading behavior in the extensional case requires 
introduction of a 'crinkle factor' describing a rapid loss of segmental ordering. We discuss in more 
detail here the physics of this, which we argue involves non-entropic contributions to the polymer 
stress, and which might lead to some important differences between shear and elongation. We also 
discuss some fundamental and possibly testable issues concerning the physical meaning of entropic 
elasticity in vitrified polymers. Finally we present new results for the startup of steady shear flow, 
addressing the possible role of transient shear banding. 

PACS numbers: 64.70.pj,62.20.-x,83.80.Va 



I. INTRODUCTION 

Understanding the flow of polymeric materials is a cen- 
tral issue in their manufacture and performance. For 
molten systems, profound insights into polymer rheology 
can be obtained by combining simple ideas on entropic 
elasticity (stress arising from entropy changes among 
random- walk chains) with conceptually simple but math- 
ematically sophisticated ideas on how molten polymer 
chains can move under strongly entangled conditions 
[U |2] ■ Despite their ubiquity in products of all kinds, 
progress in understanding the flow of polymer glasses 
has been much slower. This is perhaps not surprising 
given that even 'simple' glasses, arising in hard-sphere 
atomic or colloidal materials with no polymeric degrees 
of freedom, lead to constitutive models of roughly similar 
complexity to that of molten polymers [3J 3] . These non- 
polymeric glasses [5HZ] share with polymeric ones [51 H] 
the feature of undergoing slow plastic deformation in re- 
sponse to applied stress. This applies at least in the tem- 
perature range just below the glass transition tempera- 
ture T g which we address here. (At much lower temper- 
atures, brittle fracture may occur.) However in general 
they do not show strain-hardening as seen in polymers: 
simple glasses are melted by stress, and at constant load 
then remain indefinitely in a flowing state [10] . 

Faced with the complexity of the polymer glass prob- 
lem, a number of modeling strategies are possible. A 
crucial question is how much glass physics and how much 
polymer physics to try to include. One important avenue 
has been to first develop a largely empirical approach to 
the glass sector, coupled initially with minimalist poly- 
mers (modelled as elastic dumb-bells, for example) and 
then add more polymeric detail such as additional Rouse 



modes. Exemplary of this approach is the EGP (Eind- 
hoven Glassy Polymer) model [1202, whose antecedents 
date back to the work of Haward and Thackray [T3J . 

On the other hand, our understanding of flow in sim- 
ple glasses has improved greatly in recent years, through 
advances in microscopic [31 [TU] and mesoscopic [TH - H5] 
theory, much of it well tested experimentally. Crucial to 
glass rheology (but not yet under good control in the mi- 
croscopic approaches [3J|4]) is physical aging: a quiescent 
glass becomes more solid over time, but 'rejuvenates' un- 
der flow into a liquefied state. This physics is distilled by 
so-called 'fluidity' models, which feature a time evolution 
equation for a single structural relaxation rate, called the 
fluidity [T71[TS]. In the so-called 'simple aging' scenario, 
the inverse of the fluidity, which is the structural relax- 
ation time r, for a quiescent glass increases linearly with 
its age [9l [HI [19] . A slow steady flow cuts off this growth 
at the inverse flow rate; increasing the shear rate (e.g., 
in startup of steady shear) causes rejuvenation [T71 ITS] , 

We believe that this improved understanding of the 
glass sector justifies a renewed approach to the problem 
of polymer glasses, again starting initially by coupling a 
suitable glass model to the simplest (dumb-bell) descrip- 
tion of the polymer degrees of freedom. (Clearly, a future 
objective is to improve the polymeric part of the descrip- 
tion, but it is first instructive to see what physics can 
emerge at dumb-bell level.) This renewed approach was 
initiated by three of us in [5D], mainly within the spe- 
cific context of explaining strain hardening under tensile 
loading, and its effects on segmental relaxation rates, as 
observed in [5J. In the present paper, we give a fuller 
presentation of the same approach, and also give new 
predictions for the corresponding creep-test experiments 
in a simple shear geometry. To the best of our knowledge 



such tests have not yet been done, at least for the exper- 
imental system addressed in [8] , so this work (which uses 
only the fit parameters from the extensional data in that 
paper) offers new and separately testable predictions of 
our simple approach. We also present some additional 
results on shear startup focussing on the possible role of 
transient shear banding. 



II. STRAIN HARDENING 

In polymeric glasses, the interplay between polymeric 
and glassy degrees of freedom causes new properties to 
emerge. Particularly striking among these is the response 
to loading of the segmental relaxation time which 
controls the rate of local rearrangements. Lee et al [5Il2"T] 
found that r(t) falls steadily during the early stages of 
elongational deformation, and then falls more sharply, 
reaching a small fraction ~ 10 -3 3 of its initial level before 
dramatically rising again. (FigJT] shows all these features 
for a theoretical curve that was presented in 20] as the 
best fit to the data of Lee et al; the unloading part is 
discussed separately below.) The rise in r(t) happens as 
the local strain rate starts to drop: the latter property 
defines onset of the so-called strain-hardening regime. 

Various elements of this scenario have been confirmed 
in coarse-grained and molecular simulations 22-25], but 
we believe ours is the first simple theoretical picture of 
it to capture the key features just described. For in- 
stance the theories of Chen and Schweizer [T3 l2l)ll28| . 
which emphasize the role of stress-induced hopping over 
barriers (following Eyring [29] ). seem unable to account 
for the striking dip in segmental relaxation time. More 
generally the assumption of a stress-dependent fluidity 
sits uneasily with the idea of aging in glasses, which 
causes time dependence of the fluidity even under con- 
stant stress [HI HH HI]- In our own work, in common 
with more fundamental theories of simple glasses [3] H] , 
fluidity arises primarily by strain-induced barrier crossing 
(deformation forces cages to break) as opposed to stress- 
induced effects. This distinction means that r(i) depends 
on the accumulated flow history rather than solely on 
the present state of stress; the latter possibility is clearly 
ruled out by the data of [5] . 

Previous work to incorporate aging and flow- 
rejuvenation into polymer glass theory led to the Eind- 
hoven Glassy Polymer model [TTj . whose viscosity (or 
equivalently r(i)) is controlled by the evolution of a state 
parameter S that is age- and strain-dependent. Such a 
state parameter might have various physical meanings 
and - especially if one is also interested in the tempera- 
ture and pressure dependence which is very important for 
applications [11] - its time evolution is in general likely 
to be complicated. Given this, it is hardly a criticism of 
that work that empirical simplifying assumptions were 
made concerning the time evolution of S. 

Crucially, however, among these assumptions is one 
that cannot be reconciled with the aging and rejuve- 
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FIG. 1: (Color online) Solid curves: local strain i — expe— 1 
[8], reduced relaxation time r(t)/t m and tensile stresses T p ' s = 
G p ' s (a p ' z s — crJZ'x) of the polymer (p) and solvent (s) during 
loading of an infinite uniform cylinder. Here t w — t(0~) is 
the age of the system when the experiment begins. Param- 
eters are G s /G p = 8.5, fi = 12.5, t w /r = 10 4 , r = 6s; 
applied force / initial area / = 2.7G P . (The curve for T p , in 
red, initially lies below T" but crosses it during strain harden- 
ing.) The unload results for the basic model (9 = 1) is shown 
dashed; the solid curve after unload has 8 — 0.1. The hori- 
zontal axis is marked both in dimensionless model units (top) 
and real time (converted using To), bottom. (As explained 
in [2D], the numerical solver introduces, in lieu of inertia, 
a small additional fluid viscosity r) n = 0.05G p tq into Eq[T] 
whose magnitude has negligible influence on these plots.) 



nation scenario encapsulated by the fluidity approach. 
Specifically, the EGP model assumes aging and rejuvena- 
tion to have factorable effects on S. (The authors of [11] 
already acknowledged the potential limitations of this ap- 
proximation.) Thus, strain-induced rejuvenation causes 
a (cumulative) loss of structure, encoded in a reduction 
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FIG. 2: Schematic evolution of the relaxation time r(t) in a 
sample with r(0 - ) =t w - which is its age or "waiting time" in 
our model - subjected later to a step strain causing a sudden 
drop in r. In the simple aging picture, r(t) rebuilds from this 
point with the same slope as before (upper curve). However, 
if aging and rejuvenation are factorable, the slope of the curve 
drops by the same factor as r does (lower curve). 



in S, which is remembered indefinitely and reduces all 
subsequent relaxation times by the same factor. This 
contrasts with fluidity models, where strain-induced flu- 
idization resets the clock for aging but does not change 
its rate. Fig(2] schematically compares these two cases 
following a step strain. 

Although we have not seen any direct attempts in the 
literature to fit the data of [5] to the EGP model, it 
seems qualitatively clear from this discussion that the 
EGP model is unlikely to account for the rapid recovery 
in r(t) on entering the strain hardening regime. How- 
ever it should be possible to do this with relatively simple 
fluidity-inspired modifications to the EGP model, and in- 
deed the work begun in [20] and continued in this paper 
could be interpreted as a step in exactly that direction. 
Meanwhile though, the EGP's precepts as currently for- 
mulated are arguably unsuited to the regime of strong 
fluidization arising just below T g that was addressed ex- 
perimentally in [5j and in recent glass rheology theories 
p2 EH [HI [HI HI]- This regime may have different physics 
from the more strongly solidified regime further below 
the glass transition (though it is not clear even then that 
factorable behavior would result). The latter regime is 
of course highly relevant to in-use mechanical deforma- 
tion of glassy polymers (in car fenders and the like) but 
the near-T 9 regime is certainly also important in molding 
and other processing steps for such materials. Possibly 
both regimes could be spanned by a single, generalized 
evolution equation for S (or r) but we leave this aspect 
to future work. 



III. DUMB-BELL FLUIDITY MODEL 

Despite several recent efforts [HI [2BJ [251 130], creat- 
ing a comprehensive theory of rheological aging in poly- 



mer glasses remains a formidable task. Here we recall 
the ingredients of a minimal model |20j . combining just 
two key elements of any such theory (nonfactorable ag- 
ing/rejuvenation, and the strain dependence of polymer- 
borne stresses), which semiquantitatively explains many 
of th e results reported in [8j. Subsequently (Sections 
V|VI ) we give a more detailed discussion of the underly- 
ing physical assumptions than was possible in [20] . 

We consider polymeric dumb-bells [2] suspended in a 
glassy fluid or 'solvent' whose microscopic relaxation time 
r(t) obeys a fluidity- type equation showing simple ag- 
ing and flow-rejuvenation. The glassy fluid could be a 
genuine solvent, or alternatively, if no actual solvent is 
present, it could describe the short-scale, relatively fast 
degrees of freedom that control the local dynamics of 
segments. The model distinguishes these from the slower 
degrees of freedom of large sections of chain, represented 
here for simplicity by the added dumb-bells. This 'di- 
vision of labor' is well established in polymer melt the- 
ory where, for instance, the Rouse model adopts a local 
friction (equivalent to a time-independent segmental re- 
laxation time r) that can equally be provided by drag 
against a solvent or against segments of other chains. In- 
deed there is little difference in behavior between molten 
polymers and sufficiently concentrated polymer solutions 
[2J. One interesting question is whether this also applies 
to polymer glasses: are there important differences in be- 
havior between a single-component polymer glass, and a 
concentrated polymer solution in a molecular solvent be- 
low its T g ? Our model assumes that any such differences 
are controlled primarily by the relative magnitudes of the 
elastic moduli of polymeric and solvent-like degrees of 
freedom (see below) . Dilution of the polymer component 
would delay the onset of strain hardening and enhance 
the solvent yield behavior that precedes this. However, 
we do not pursue this issue in the current paper. 

For simplicity we treat the dumb-bells first as linearly 
elastic elements - as is valid for the entropic elasticity 
arising at modest deformations in the molten state [2J. 
However, in Sections |V|VI wc discuss further the true 
nature of the polymer stress in polymeric glasses which 
is not solely entropic in character j28j [36] , and show how 
this can be partially incorporated into the model. 

We start by defining a deviatoric stress tensor 



G P ((T P 



G s 



(1) 



where <t p and cr s are dimensionless conformation tensors 
for polymer and solvent, G p,s associated elastic mod- 
uli (see below), and I the unit tensor. For entropic 
dumb-bells of concentration n and equilibrium mean 
square end-to-end distance (R 2 )o, G p = 3nksT and 
<j p = 3(RR)/(i? 2 )o. The dumb-bell conformation tensor 
o p is taken to obey a standard upper-convected Maxwell 
model, which can be derived [2] by considering the dy- 
namics of an ensemble of dumb-bells whose endpoints are 
(i) advected by flow (ii) subject to a linear spring force 
inwards along the line connecting them and (iii) subject 
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to independent Brownian diffusion: 

& p + v.VcrP = <x p .Vv + (Vv) T .crP - a{(T p - I)/r (2) 

Here the structural relaxation time t p — r/a is propor- 
tional to, but much larger than, that of the solvent, which 
we denote r. In the simplest models of dense, molten, but 
unentangled polymers, a = N~ 2 with N the polymeriza- 
tion index [5] . On the other hand in a lightly crosslinked 
elastomeric network, as studied experimentally in [S] , one 
expects a = 0. Note that so long as a is small enough, 
its precise value has very little effect on our numerical 
calculations, at least for the experimental protocol of [5]. 



This feature is explained in Section IV 



We now turn to the solvent stress G s cr s . Bearing in 
mind that (unless a true solvent is present) this stress 
is also polymeric in origin, we take the relevant confor- 
mation tensor to obey another upper-convected Maxwell 
equation, now with relaxation time r: 

& s + v.Ver 5 = <t s .Vv + (Vv) T er s - (a s - I)/r (3) 

Because local degrees of freedom outnumber the chain- 
scale ones we expect G s > G p . 

In a departure from molten polymer models, in which 
the Maxwell time is a fixed quantity, we next assume 
this structural relaxation time to have its own dynamics, 
governed by a fluidity equation 



f + v.Vr = 1 - (r - t )A 



(4) 



Here we have defined the following scalar invariant mea- 
sure of flow rate 



A(D) = /V2Tr(D.D) 



(5) 



with v the fluid velocity and D = (Vv + (Vv) T )/2. 
The parameter /i is an order unity dimensionless quan- 
tity which controls how effectively flow prevents aging. 
This in turn sets the steady-state relaxation time as a 
function of strain rate via 



TO 



(6) 



Without flow, r increases linearly in time at a solidifi- 
cation rate f(D = 0). This quantity is dimensionless and 
for simplicity we set it to unity. An alternative (given 
that direct measurements of the aging of segmental re- 
laxation rates were not made in [5]) would be to leave 
it in the model as a freely floating parameter. We have 
refrained from doing this, but would expect other choices 
(replacing the 1 in Eqj4]with some other number) to give 
quantitative rather than qualitative changes. 

With flow present and aging hypothetically switched 
off, t would undergo deformation-induced relaxation to- 
wards To which is a 'fully rejuvenated' (microscopic) 
value. According to Eqj3] this relaxation occurs at the 
rate A. In steady shear, A = and r varies inversely 
with strain rate 7 in accord with microscopic theory 



[3J[3T]. For uniaxial elongation at strain rate i, ^ re- 
duces to A = //V3|e|- Note that in this simple fluidity 
model, the rejuvenation of r is essentially strain-induced 
[TO] but, in contrast to the factorable model of [UJ, can 
be rapidly reversed by subsequent aging. 

Our chosen fluidity model can in principle be embel- 
lished with additional parameters, so as to incorporate 
further relevant information about glass rheology. We do 
not pursue this in detail here, but mention one interesting 
example, which is to write instead of Eq[4] 



f + v.Vr = I — rA(r — tq) — 



(1 - r)r 
A^+to 



(7) 



For our chosen model, r = 1, the dynamics of r is essen- 
tially strain-induced, so that a pulse of fast flow (whose 
limit defines a 'step strain') is just as effective at flu- 
idizing the system as when the same strain is applied 
more slowly. In contrast, for r — 0, step strains have 
no fluidizing influence at all. (This model of fluidity was 
explored in [32] ■ ) Both limits give the same result, Eqj6j 
for steady-state flows. One expects the true behavior of 
glasses to lie between these two limits — but in a manner 
that can itself depend on the strain amplitude, and sign, 
in relation to the previous deformation history [3J [33] . 

To complete the specification of our model, we add 
the standard equations of mass and force balance for an 
incompressible fluid of negligible inertia: 



V.v = 
V.[S + 2?7D] = 



(8) 
(9) 



In ([9]) we have added a small additional Newtonian vis- 
cosity r\. This is included solely for numerical reasons: 
it ensures that a state of force balance can always be 
found without fluid acceleration, so that inertia can be 
set to zero from the outset. The alternative is to intro- 
duce a small inertia; in practice for the flows considered 
here, this choice (along with the precise value of rf) has 
no discernible effect on any of the results. 



IV. ELONGATION AL LOADING 

In |20j we presented numerical data for uniaxial ex- 
tension flows within a lubrication approximation appro- 
priate to long cylindrical samples undergoing homoge- 
neous deformation. As demonstrated there, the data re- 
plotted again here in Fig. [TJ for both the 'local strain' 
e{t) = expe — 1 and the segmental relaxation time r(t) 
closely resembles the experimental results of 8J for a 
cylindrical sample of polymer glass subject to a sudden 
tensile loading introduced after a period of aging at rest, 
up to the point where the load is removed (but not be- 
yond, see Section |v| below) . 

Because the sample was lightly cross-linked, a was set 
negligibly small; the experimental protocol of [5] then 
specifies the engineering stress / = 16MPa (and also the 
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time t u — 9400s at which unload later occurs) and mea- 
sures directly the initial relaxation time r(0) = 6 x 10 4 s. 
There remain four material parameters in our model: 
G p ,G s ,t and p,. Of these, G p ~ 6MPa was deduced 
from the deformation in the strain-hardened regime just 
before unload; G s ~ 50MPa and fj, ~ 12.5 were in turn 
estimated from the step-change in r during initial load- 
ing, and from the separately measured slope [5] of the 
parametric 'effective flow curve' e(r). Hence the only 
unconstrained parameter in fitting the dip in r(i) was 
To, which was found to obey t ~ 6s. 

Striking features of the experimental data, reproduced 
by our simple model, include: the tenfold initial drop 
in t on applying the load; its subsequent further decline 
to a sharp minimum T m - ln ~ 10~ 3 ' 3 r(0) near the point 
of maximum elongation rate; and its rapid rise from that 
minimum towards a strain-hardened plateau. Other com- 
parisons using the same model parameters are given in 
[20) and also offer encouraging agreement between model 
and experiment. 

Fig. [I] shows (alongside strain and relaxation time) the 
tensile stresses T p,s carried by polymer and solvent re- 
spectively. Quite striking is the transfer of load from 
solvent to polymer which ushers in the strain harden- 
ing regime. Indeed, according to our simple model, it 
is this transfer that explains strain hardening. Within 
the model, what matters is that the glassy solvent has a 
certain yield stress Ey, initially exceeded by the applied 
load. After an initial step-down in r caused by the step 
strain on loading, the material starts to flow continuously 
which fluidizes it further. As a result its strain rate ac- 
celerates, giving positive feedback so that r(t) plunges 
downwards. As deformation builds up, however, an ever 
growing share of the applied stress is carried not by the 
solvent but by the stretching polymer chains, which have 
little or no relaxation on this timescale. This polymeric 
stress causes the flow rate to drop: the solvent, whose 
stress now obeys E s < Ey, therefore starts to solidify. It 
enters a simple aging regime in which the flow progres- 
sively grinds to a halt. 

It is notable that this explanation (which also directly 
explains the remarkable behavior seen for r(t)) does not 
require any nonlinearity in the polymeric response. (Such 
nonlinearity might nonetheless be present; see Section M 
below.) It does require the existence of two parallel chan- 
nels for supporting stress, so that the unloading of the 
solvent allows its solidification at late times. Accord- 
ingly, we predict that strain hardening of this character 
will never emerge in a simple, one component (molecular 
or colloidal) glass. 

An interesting consequence of this scenario is that the 
polymers, whose Maxwell time scales with that of the 
solvent but is much larger, can never reach a state of ter- 
minal relaxation once strain hardening has begun. Even 
if (departing from the protocol of [8]) the load were 
never removed, once the polymers are carrying most of 
the stress, the solvent re-freezes and, in accord with the 
simple aging picture, the polymer relaxation time starts 



growing as t p ~ t/a. Since a is large, at no stage can 
the terminal regime (i > t p ) be reached. This is why, as 
mentioned in Section |III[ the value of a is almost irrele- 
vant to our predictions, so long as it is small enough. Be- 
cause the terminal zone is not reached, this experimental 
protocol will not distinguish a cross linked network from 
unlinked chains of high molecular weight. 

As illustrated in Fig.[2j models that factorize aging and 
rejuvenation effects are seriously challenged by the rapid 
recovery of r after the dip. (A multimode spectrum |12j 
is unlikely to help here.) With simple aging, such factor- 
ization predicts r ~ (t+t(0))/(e), so that if the segmental 
relaxation times falls from its pre-deformation value i(0) 
to a small value r = ft(0) — r m i n at the dip, a tenfold 
recovery to r ~ 10r mi „ does not occur until t ~ 10i(0) 
which for the experiment of [8] means t~6x 10 5 s. This 
prediction is at least two orders of magnitude longer than 
the experimental timescale (~ 3000s), and is an under- 
estimate since the dip occurs during a period of rapid 
flow (so aging is less fast than in a quiescent system) . In 
contrast, for our simple- aging fluidity model, the recov- 
ery time is predicted to be of order t ~ 10/t(0) = 300s. 
This is also an underestimate, for the same reason, and 
on that basis is in better accord with the data. 



V. ELONGATIONAL UNLOADING 

A significant shortcoming of the model becomes appar- 
ent when the sample is unloaded. Here the experiments 
show a modest drop in r immediately on removing the 
load, followed by a gradual recovery towards the pre- 
deformation value. The dotted line in Fig. [T] shows the 
prediction based on Eqs. ([2}j5]); r drops, but then falls 
much further before recovering. (The solid line is a mod- 
ified model, discussed below, which much more closely 
resembles the experimental data [2D].) Moreover, for 
stresses well above the yield stress, this re-solidification 
only occurs after an almost complete recovery in the 
strain. 

This means that the model fails to capture a major as- 
pect of polymer glass behavior, which is that the glass can 
be deformed plastically and will then hold its shape with 
only modest relaxation once the load is removed. The 
problem arises because in the strain-hardened regime, 
the polymers carry a large elastic tensile stress, which 
generally exceeds Ey. 

Upon unloading, this polymeric stress is unbalanced 
by the external traction and therefore acts backwards on 
the vitrified solvent, causing it to yield. The resulting 
evolution of r(i) resembles a repeat of the initial loading 
experiment, but run in reverse; only when the polymer 
stress (now equal and opposite to the solvent stress) falls 
below Ey does revitrification of the solvent occur. By 
this time, for typical parameters, the sample has almost 
recovered its original shape. 

One possible candidate for the discrepancy with ex- 
periment lies in the response of r(t) to imposition of step 
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strains, as arise at both the onset and the removal of 
the extensional load. As described following EqjTJ this 
can be subtle for glasses (see [53]) and one could imag- 
ine circumstances in which curtailing the sharp drop in r 
immediately on unloading would have a strong effect on 
what happens subsequently. For this reason we have nu- 
merically explored the effect of reducing r in EqJTJon the 
post-unload behavior. For the basic model (r = 1) the 
post-unload drop in r(t) is by a factor of order 20 (for the 
model parameters used previously, see Fig.JT]); this factor 
can indeed be reduced by dropping r (the experimental 
reduction factor is of order 2) in accord with the discus- 
sion of step strain response following Eq[7j However, this 
adjustment cannot prevent the near-complete strain re- 
versal following unload, and the argument in the previous 
paragraph (which does not depend on r) explains why. 
The incomplete strain reversal can only be explained if 
the polymer stress falls below the solvent yield stress af- 
ter only a modest degree of reversal has taken place. 

Thus we are led to seek a mechanism that can allow 
the polymer stress to drop much faster on unloading than 
would normally be the case for elastic dumb-bells. An im- 
portant clue is that the value of G p ~ 6 MPa needed to 
fit the loading data is one order of magnitude larger than 
the rubbery modulus of the same material above its glass 
transition (see, e.g., [35]). This accords with a widely 
held literature view that the strain-hardened modulus 
of polymer glasses does not primarily arise from single- 
chain entropic elasticity [28J [36]. The larger modulus 
could come from intrinsic nonlinearity (finite extensibil- 
ity) in the stress-strain curve of a single chain or it could 
be caused by the buildup of local mechanical interaction 
forces at chain-chain contacts. These two mechanisms 
are much more similar than they might first appear. In- 
deed, for an inextensible rod in an elongational flow field, 
the stress is primarily a viscous solvent stress and caused 
by the obligatory distortion of fluid streamlines to main- 
tain a no-slip boundary condition with the rod. Given 
that we are treating the local friction between chains by 
pretending each chain is embedded in a glassy solvent, 
the chain-chain contact forces can likewise be thought of 
arising within the solvent. 

As proposed in [20] . our currently preferred explana- 
tion invokes a well established piece of physics that arises 
when flexible polymer chains are placed in an extensional 
flow fast enough for them to extend rapidly relative to 
their own relaxation time. This requires et p 3> 1; within 
our model and for well aged samples, this condition al- 
ways holds when a is small. For instance, in steady state 
t p ~ fii/a so that the fast flow regime applies, indepen- 
dent of the actual strain rate, whenever a <C n ~ 0.1. 

In that regime it is argued that short subsections of 
chain quickly stretch close to full extension locally, form- 
ing a quasi-one-dimensional filament containing hairpin- 
like kinks [3JJ [38] . Further stretching is mediated by mi- 
gration and annihilation of neighboring kinks of opposite 
sign. During this process, much of the stress carried by 
the polymers is not entropic-elastic, but instead caused 



by viscous drag against the extended subsections of chain 
(which behave locally as inextensible rods). 

If one tries to describe this situation with a dumb-bell 
model, then, to match the resulting stress, an enhanced 
value of G p must be invoked. However, if the flow is now 
suddenly stopped (or reversed) much of this extra stress 
disappears on a relatively rapid timescale [37] . This is not 
set by the global relaxation time of the whole chain, as 
a dumb-bell picture would assume, but the time it takes 
for locally straight subsections to crinkle back into an 
entropic conformation. A similar physical process arises 
when chains are fully stretched at all scales, again causing 
a rapid loss of polymer stress on unloading with only 
partial relaxation of global polymer conformations 39 . 

A proper treatment of these effects requires explicit 
analysis of short-scale polymeric degrees of freedom, cor- 
responding to higher Rouse modes rather than the sin- 
gle mode retained in the dumb-bell description (3[ [37] . 
A phenomenological 'renormalized dumb-bell' represen- 
tation can however be gained by assuming that, on un- 
loading, the effective polymer modulus abruptly drops 
by a 'crinkle factor' 9 so that G p -» 9G P . Here the 
lower value should roughly correspond to the standard 
entropic-elastic modulus [3D] . 

A numerical value 9 = 0.1 for the crinkle factor is sug- 
gested by the previous observation that G p is about ten 
times larger than the entropic modulus observed above 
T g . The solid lines in Fig. [I] show the results of this 
choice, which are much closer to the experimental be- 
havior both for the strain and for r(t). Post-unload, the 
polymer stress acting backwards on the solvent is now 
safely below the solvent yield stress; the result is a mod- 
est drop and then slow increase in r(t). This is accompa- 
nied by only partial restoration of the initial shape with 
arrest in a finite state of deformation when the entropic 
polymer stress 9G p cr p falls below the yield threshold for 
the glassy solvent. 



VI. SOME ISSUES OF PRINCIPLE 

Given the crudeness of the crinkle-factor picture, this 
appears a quite satisfactory outcome. The approach does 
highlight certain questions of principle however. For in- 
stance, it is sometimes argued that entropic contribu- 
tions to thermodynamic quantities such as stresses can 
only arise under conditions where kinetic exploration of 
configuration space actually occurs [41] . Holders of this 
view might expect the entropic part of the polymer stress 
to switch off as the solvent freezes, requiring an extra 
creep of the sample during freezing until the non-entropic 
polymer stress has increased enough to replace the miss- 
ing entropic contribution. We argue instead that the 
trapped polymer conformations continue, microstate by 
microstate, to exert exactly the same mechanical force as 
they did when exploring the equilibrium ensemble. Since 
the thermodynamic tension in a chain is also the time- 
average of a thermally fluctuating mechanical tension, 
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a 'purely' entropic polymer stress can indeed be frozen 
in on vitrification. (However, if the sample is deformed 
again subsequently, this part of the stress has no reason 
to evolve as predicted by entropic elasticity, since poly- 
mer conformations are trapped with the 'wrong' proba- 
bility distribution for the new conditions.) 

A viewpoint that frozen polymers cannot exert an en- 
tropic stress would render our strain-hardening mecha- 
nism highly problematic in any system where the crinkle 
factor mechanism was not operative, so that the poly- 
mer stress is purely entropic throughout. (This might 
arguably include the case of shear deformation, consid- 
ered next.) For in the strain-hardened state as t — > oo, it 
is only because the polymers carry the stress that the sol- 
vent can remain frozen. Clearly if solvent freezing causes 
the polymer stress itself to disappear, then we have a 
paradox, and something else has to happen, such as a 
solvent stress that asymptotically approaches £ y from 
above as t — ¥ oo. The paradox is avoided if we are cor- 
rect that an entropic polymer stress can persist within 
a frozen solvent. Note that closely analagous arguments 
concerning the persistence of entropic contributions to 
the chemical potential in frozen systems have been con- 
firmed in both experiment |42] and simulation [43] . 



VII. SHEAR DEFORMATION: STEP STRESS 

Experimentally a tensile loading test as presented in 
[5] is a standard procedure for polymer glasses. The ge- 
ometry is quite simple initially (a cylinder clamped at 
both ends) but the fact that the cross-section narrows 
with time means that the actual stress deviates from the 
engineering stress in the material and is not constant in 
time. This situation is made worse by 'necking' which is 
commonplace in such tests and was ignored in the dis- 
cussion above (but is briefly addressed within our model 
in [20]). 

From a theory or modeling perspective, it is interesting 
to ask what happens when a similar application and then 
removal of a constant load is made in a simple shear 
geometry. From the physics that was described in Section 
IV based on the precepts of standard dumb-bell elasticity, 
one expects in principle similar type of strain hardening 
to arise. The melting of the solvent at shear stresses 
above its yield stress, the gradual transfer of stress from 
solvent to polymer, and the subsequent revitrification of 
the solvent should all happen in a similar way. 

We therefore plot in Fig.|3]the predicted response from 
our model for a shear load and unload protocol. This uses 
the same material parameters as in the fit to the experi- 
ments of [8] and is at a matched value of the engineering 
stress and the unload time to allow direct comparison 
with the elongational case as presented in Fig. [T] As 
anticipated above we find broadly similar behavior al- 
though the stress dip is far less pronounced. This latter 
effect may be due to the reduced strength of feedback 
during the initial shear-melting phase: recall that an ad- 
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FIG. 3: (Color online) Solid curves: shear strain 7, reduced 
relaxation time r(t)/t m and shear stresses E p,s = G p ' s a^y of 
the polymer (p) and solvent (s) during shear loading. Pa- 
rameters as in Fig[T] with a matched ratio of shear stress to 
solvent yield stress. (The curve for E p , in red, initially lies 
below S s but crosses it during strain hardening.) The unload 
results for the basic model (9 = 1) is shown dashed; the solid 
curve after unload has 6 = 0.1. 



ditional geometric feedback (strain increasing the stress 
due to cross section change) is present only in elongation. 
Perhaps for similar reasons the dip, which is again syn- 
chronous with the stress transfer from solvent to polymer, 
appears somewhat later in shear than in elongation. In 
the baseline model with 9 = 1 (as might describe purely 
entropic elasticity) the second dip in r(t) immediately af- 
ter unloading is also less pronounced than in elongation. 
Introducing a crinkle factor 9 = 0.1 has similar effects to 
what happens in elongation, making this second dip still 
less pronounced and, perhaps more importantly, drasti- 
cally reducing the strain recovery achieved after unload- 
ing. Fig. [4] shows the corresponding data for a substan- 
tially higher stress (six times the yield value). This is the 
value needed for the initial dip to match that seen in elon- 
gation (a factor 10~ 3 ' 3 in t/t(0)); unsurprisingly it now 
happens very much earlier since the deformation builds 
far more quickly at this high stress. In this run we also 
reduced the unload time t u to roughly match t(£ u )/t(0) 
at the point of unload. 

To ensure comparability with the elongational results 
presented in Fig. [T] and [20] , the results of Figs 3j4 have 
been calculated under conditions where homogeneity of 
the flow is imposed. (The numerical protocols are also 
as detailed in 20 .) An interesting question is whether 
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FIG. 4: (Color online) As in Fig. [3] but with shear stress 
increased by a factor 2.76. (This roughly matches the depth 
of the minimum in r(t) to the elongational data of Fig. [I]) 



the step stress loading protocol can create strain inhomo- 
geneities such as shear bands (layers of large and small 
shear strain coexisting at equal stress). A theoretical 
criterion for shear banding in step stress has recently 
been derived; this requires that the underlying homoge- 
neous creep curve j(t), differentiated to give j(t), obeys 
dfi/dti > [33]. We have explored this question nu- 
merically using the same methods as explained in Section 
|VIII| below for shear startup. We found that including 
noise in the initial condition can create persistent strain 
inhomogeneities throughout the loading and unloading 
phase. However, although these do not decay mono- 
tonically in time, they remain extremely small through- 
out. For example, with parameter values chosen as in 
Fig. [4j a fractional variation in strain rate, initially of or- 
der 10 -3 , decayed rapidly before rising slowly again to a 
peak whose height was itself less than 10 -3 . (The height 
of this peak was found to be linear in the initial noise 
amplitude.) The corresponding strain inhomogeneities 
remain smaller still for the more modest stress applied 
in Fig. [3j Moreover, we found that increasing the mod- 
ulus, by curtailing the growth in shear rate through the 
strain hardening effect, suppresses strain inhomogeneity 
in step-stress relative to cases where G p is small but fi- 
nite. 

All this contrasts somewhat with the case of elonga- 
tion where necking (ignored above but addressed in [3U]) 
is a ubiquitous phenomenon in polymer glasses. It stems 



mainly from the feedback between cross section and local 
stress (so that a region of narrowing in a cylindrical sam- 
ple then narrows further). It also contrasts with expec- 
tations based on recent work in startup of steady shear 
(see Section VIII ) , which does give shear banding, albeit 
transient, in fluidity models without polymer |32j . 



Before addressing the case of startup (Section VIII ) we 
return to the crinkle factor 0. The reasoning for 9 < 1 
given in Section [V] was based on the picture of locally 
stretched chain segments creating high viscous dissipa- 
tion |37[ 138] ; this picture was originally developed in the 
specific context of elongational flows. Shear flow of course 
also has an elongational axis, but for large strains this is 
not aligned with the polymer orientation (the latter ap- 
proaches the flow direction while the elongational axis 
remains at 45 degrees). For this reason, within the pic- 
ture we have developed, one can perhaps expect a much 
weaker buildup of non-entropic stresses and hence a much 
smaller drop in polymer modulus on unload. Indeed, 
taken at face value our arguments suggest a value of 9 
close to unity, which in turn requires that the polymer 
modulus relevant to shear flow is the entropic one, and 
hence approximately ten times smaller than that used to 
fit the elongational loading data [35]. 

If this grossly oversimplified picture were correct then 
one could expect very significant differences in the exper- 
imental curves, not only for r(t) but also for strain, be- 
tween elongational and shear loading. Crudely speaking 
(i.e., ignoring the required shift in G p ) the elongational 
experiment would follow the solid line in Fig. [T] as indeed 
it nearly does, but the shear protocol would follow the 
dotted line in Fig. [3] The most important macroscopic 
consequence is that in shear the strain recovery on un- 
load would be much larger, because the polymer loading 
on the solvent, which is nearly the same in magnitude as 
the initially applied load, is always enough to re-melt the 
solvent (unless the initial load was only slightly above the 
solvent yield stress). 

To counter this argument, it has recently been sug- 
gested that strain hardening in polymer glasses is at- 
tributable not to the viscous stress on a kinked fila- 
ment but nonetheless to the 'pulling tight' of a small 
fraction of chain strands trapped between entanglements 
[46] . This is quite similar in spirit to the mechanism dis- 
cussed above, and could also lead to small 9, but could 
be somewhat less dependent on the deformation geom- 
etry. To help understand these issues further, it would 
be extremely interesting to see experimental comparisons 
between tensile and shear loading on similarly prepared 
samples at comparable stress. 



VIII. SHEAR STARTUP 

To check that our simple model is not limited to the 
step-stress case but also behaves reasonably in strain- 
controlled flows, we calculated in [20J (as supplemental 
material) the stress responses for startup of steady elon- 



9 



gation and compression. These show a stress overshoot, 
whose height varies as hi(it w ), which is similar to the 
behavior found (under shear) in simple aging fluids [13] , 
and also broadly accords with reports in the polymer 
glass literature (see e.g. [TT1IT2] . 

Here we extend this further to address the startup of 
steady shear which is a widespread experimental pro- 
tocol in the polymer rheology literature generally. The 
physics of the stress overshoot involves a transient elas- 
tic deformation that exceeds the recoverable strain in the 
steady state that is reached eventually. This is common 
in systems with strong shear thinning whose initial elastic 
response involves structures, or polymer conformations, 
that can be reached at modest strains from the initial 
state but cannot be sustained thereafter. 

Figure [5] shows the polymer, solvent and total contri- 
butions to the shear stress during startup for parameters 
as in Fig. [3] at an imposed strain rate 7 = 10~ 2 . These 
curves are calculated numerically on the assumption that 
the flow remains homogeneous. 

However, recently it has become clear that overshoots 
of this kind, which are particularly common in aging sys- 
tems, often cause instabilities leading to transient shear 
banding (layers of material with different 7 at a com- 
mon stress) [37] • These instabilities typically 
occur within the region where stress is decreasing and 
strain increasing, which can crudely be viewed a nega- 
tive differential elastic constant (although this picture is 
certainly oversimplified 44 ). We therefore repeated the 
simulations allowing the flow to become nonuniform and 
indeed find transient strain rate variations with position 
that rise to a maximum near the stress peak and then 
decay only very slowly; these results are also in Fig. [5] as 
is the total stress arising when the nonuniformity is taken 
into account. (The numerical methods are summarized 
in @S].) 

We should emphasize that for these parameters (chosen 
to match the elongation experiments of jS]), the shear 
banding is a very weak effect. Although the system is 
transiently unstable to the formation of inhomogeneities 
in the velocity gradient, these barely deserve to be called 
shear bands. Indeed the maximum deviation in shear 
rate plotted in Fig. [5] depends directly on the initial noise 
used to excite the perturbation. For the data shown the 
velocity inhomogeneity in the initial state (which decays 
rapidly in the first few timesteps) was in fact larger than 
the maximum recorded subsequently. 

However, the banding does become robust and repro- 
ducible at much larger values of t w /ro where we recall 
that t w = t(0~) is the sample age at the initiation of 
shear. For one such case (t w /ro — 10 s ), this is demon- 
strated in Fig. [6] which also shows a comparison between 
our polymer model and a simple fluidity model obtained 
by setting G p to zero. A linear stability analysis [35] 
shows that the flow becomes unstable, as expected, in 
the decreasing part of the overshoot but that the ad- 
ditional polymer stress, which has no overshoot, exerts 
a slight stabilizing role. To confirm this we studied the 
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FIG. 5: (Color online) Startup of steady shear at applied 
strain rate 7 = 10~ 2 for a system with the parameters of 
Fig. [I] Upper panel: Total shear stress (solid) and the poly- 
mer (dotted) and solvent (dashed) contributions when the 
flow is imposed to be uniform. The additional dotted points 
barely distinguishable from the solid curve show the total 
stress after shear banding is allowed for. Middle panel: the 
'degree of banding' (found by subtracting the smallest from 
the largest shear rate present at any time) for a 2D run with 
heterogeneity allowed. (A small diffusivity was added to the 
governing equations for all stress components and for r, and 
the system was initialized with a small spatially varying noise; 
see [IE]-) Lower panel: snapshots of the strain rate as a func- 
tion of position y in the flow gradient direction with symbols 
identifying strain values as in the middle panel. 



real part of the largest eigenvalue in the stability analysis 
as well as the resulting strain variation amplitude; both 
decrease monotonically with G p /G s . Despite this, the 
presence of polymer does add significant complication to 
the deviations from uniformity that actually arise, whose 
magnitude now shows an oscillation in time that was not 
present in the simple fluidity model (Fig. |6|. 

In shear startup, the general effect of increasing the 
polymer modulus is to decrease the degree of transient 
banding observed. This is because the polymeric stress 
contribution increases monotonically with strain (there is 
negligible polymer relaxation) which decreases the sever- 
ity of the stress decline after the overshoot in the total 
stress. Because transient shear banding is triggered by 
the declining stress |33j , it is mitigated by the polymeric 
contribution. 
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FIG. 6: (Color online) Upper panel: shear stress as a function 
of time in shear startup for the nudity model without (lower 
curve) and with polymer. (G p /G s = 0,1/8.5 respectively.) 
Dashed regions: transient instability as found by linear sta- 
bility analysis |48j . Dotted curve is total stress allowing for 
inhomogeneous flow. Middle panel: the resulting absolute 
value of shear rate variations with polymer (multiply bumped 
curve with added symbols) and without (singly cusped curve) 
[48| . Lower panel: snapshots of strain rate profile with poly- 
mer present, at strain points identified by the symbols as in 
middle panel. Note the strong strain inhomogeneity (shear 
banding) at strains just beyond the stress overshoot. Param- 
eter values G S ,T = 1, ju = 12.5, t w = 10 s , 7 = 10" 3 . Other 
parameters as in Fig. [T] and [3S] . 



IX. CONCLUSION 

In this paper we have developed a simple approach, ini- 
tiated in [50], to the modeling of polymer glasses. This 
builds on recent models of rheological aging and rejuvena- 
tion in simple glassy fluids, coupled to a minimal model 
of polymers at dumb-bell level. The presence of single 
polymer mode rather than a Rouse spectrum means that 
we cannot address the non-exponential form of local re- 
laxations that was reported in [5]; however it would be 
simple in principle, if cumbersome numerically, to add 
additional modes to our picture. An alternative route to 
nonexponential relaxation is through the solvent sector 
where the fluidity model similarly replaces a complicated 
relaxation lineshape (ubiquitous in simple glasses close 



to T g ) with a single relaxation time. Improvements are 
probably achievable without leaving the fluidity frame- 
work but only at the cost of additional parameter fitting. 

That issue apart, our model can explain much of what 
happens experimentally [8! when a polymer glass is sub- 
jected to elongational load. The unloading behavior is 
more problematic, but consistent with a plausible modi- 
fication of the same model which invokes on unload a re- 
duction in polymer modulus via a 'crinkle factor' 9 < 1. 
This crudely allows for the presence of non-elastic poly- 
mer stresses that are known to arise whenever et p is large 
We tentatively associate these stresses with the 
observed uplift in the polymer modulus G p beyond the 
entropic value G e that would arise for non-glassy poly- 
mers, so that 8 ~ G e /G p . However, this relies on physi- 
cal reasoning that was so far developed in the literature 
primarily to address elongational deformation, not shear. 

Partly with this in mind, and in the hope of stimulating 
new experiments, we have presented here analogous cal- 
culations for a shear loading geometry. On loading, these 
predict very similar behavior for the strain and the seg- 
mental relaxation time r(t) although the characteristic 
'dip' in the latter quantity is somewhat less pronounced. 
The unloading behavior is again similar if the same crin- 
kle factor is used. However, it is at least conceivable that 
9 in this non-elongational case should be close to unity. 
If so, the model predicts a very different unloading re- 
sponse in which the polymer stress is enough to remelt 
the solvent, causing a drastically larger strain recovery 
than seen in elongation. Such an interpretation would 
also imply a much smaller non-entropic uplift in polymer 
modulus for shear flow than for elongation. 

We have also presented new results for startup shear 
flows and shown that, in common with other shear- 
thinning glassy materials, transient flow inhomogeneities 
can be expected in this case. However these mainly stem 
from the glass physics of the solvent itself, and in fact 
the polymers have a mild stabilizing effect on such in- 
stabilities. In the polymer glass context they may there- 
fore prove inessential, although their possible presence 
should certainly be borne in mind when interpreting ex- 
periments. 

Our work suggests that an accurate representation of 
aging and rejuvenation physics may form a key part of 
any more comprehensive theory of polymer glass rheol- 
ogy. It also suggests that insights gleaned from studying 
the rheology of simple glasses will form a useful input 
into such a theory. As shown in this paper, some of these 
insights (such as simple aging and shear rejuvenation) 
can be embodied in models that are simple enough to be 
coupled with standard, if simplified, polymer rheology 
models such the dumb-bell representation. In general it 
seems a good strategy for future work to improve the 
glass and polymer sectors in tandem, rather than let one 
become much more sophisticated than the other. More 
generally our hope is that a more comprehensive account 
of polymer glasses can be achieved by judiciously com- 
bining existing types of nonlinear theory, describing non- 
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glassy polymers and simple glasses respectively. 

Of course, it may be that polymer glasses involve 'new 
physics' that is not present in either such description, 
but it may also be that many of their properties emerge 
simply from the combination. Strain hardening, accord- 
ing to our interpretation, is a case in point. The task of 
combining polymer and glass theories in this way is not 
necessarily a simple one and could raise various issues of 
principle. We have discussed one of these in this paper, 
namely the question of what is meant by entropic elas- 
ticity in a system whose conformational state is frozen. 
We believe we have the right answer to this question (the 
entropic force is maintained but its response to further 
deformation is altered), which mirrors similar discussions 
of phase equilibria involving nonergodic materials that 
date back over 80 years [321 ESI EH] ■ This type of ques- 
tion might remain a technical nicety if the predominant 
polymer stress in glasses is never entropic [30l ES] , but 



comes to the fore if the opposite is ever true. 
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